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Abstract 

In many models of Gamma-Ray Bursts (GRBs) relativistic shells are responsible for the overall envelope 
of emission. We use kinematics and symmetry to calculate the time history and spectral evolution expected 
from a relativistic shell including effects from intrinsic variations in the shell's intensity and spectra. We 
find that the decay phase of an envelope is produced by photons delayed by the shell's curvature. These 
delayed photons are produced by regions that are off-axis such that the spectra evolve according to a 
universal function (oc T~^) regardless of intrinsic variations in the rest frame of the shell. We compare 
our predictions to the overall envelope of emission of GRBs. The observed spectra evolve faster. Intrinsic 
variations cannot make the spectra evolve that fast, which adds strength to the "shell symmetry" problem: 
models, in particular, the external shock model, that involve relativistic shells must either confine the 
material to narrow pencil beams, be very inefficient, or break the local spherical symmetry so that the 
shell acts like a parallel slab. In the case of the internal shock models involving winds (i.e., central 
engines), it will probably be easier to break the local spherical symmetry, but then one must postulate 
nearly continuous energy generation at 10^^ erg s^^ lasting up to hundreds of seconds at the central site. 
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1. Introduction 

Gamma-ray Bursts (GRBs) have defied explanation 
since their discovery 25 years ago (Klebesadel, Strong, 
& Olsen 1973). The mystery was deepened by the ob- 
servation by the Burst and Transient Source Experiment 
(BATSE) that GRBs are isotropic, yet there are fewer 
dim bursts than expected relative to the bright bursts 
(Meegan et al. 1992). Two conflicting explanations 
have arisen. In the "cosmological" explanation (see, e.g., 
Paczyhski 1995), bursts originate at redshifts the order of 
1 and literally flood the entire universe with gamma- 
rays. In the alternative "galactic" explanation (see e.g., 
Lamb 1995), the bursts originate on high- velocity neu- 
tron stars which have traveled sufficient distances (~ 100 
kpc) from the disk of the Galaxy that we appear to be 
at the center of their spatial distribution. 

In both explanations, the lack of apparent photon- 
photon attenuation of high energy photons ('^ 1 Gev) 
implies substantial bulk relativistic motion. The bulk 
Lorentz factor, F, must be the on order of 10'^ to 10'^ 
for the cosmological explanation, and 10 to lO'^ for the 
galactic explanation. 

Two major scenarios involving relativistic shells have 
been developed. In the "external" shock models (see, 
Meszaros & Rees 1993, Piran, Shemi, & Narayan 1993), 
a relativistic shell is generated by a single, short-duration 
explosion at a central site. The shell coasts in a gamma- 



ray quiet phase for a time to- Eventually, the shell be- 
comes gamma-ray active, perhaps due to shocks formed 
as the shell sweeps up the ISM. (Hence, the name "exter- 
nal" shock.) Variations in the shell explain the subpeaks 
observed in GRBs, but the overall envelope of emission 
should follow that for a single expanding shell. There is 
a single, very thin, relativistic shell, and its properties 
predict the overall envelope of emission of the GRB. 

Instead of an impulsive, single injection of energy, "in- 
ternal" shock models (Rees & Meszaros 1994) feature a 
"central engine" that releases energy on timescales com- 
parable to the observed duration of the GRBs (up to 
hundreds of seconds). In some models, the central en- 
gine produces a relativistic wind which develops insta- 
bilities that generate gamma-ray producing shocks. In 
others, it produces a series of relativistic shells, and the 
fast shells collide with the slow shells resulting in a shock. 
In either case, the emitting region need not be a single, 
thin shell as in the external models. The duration of the 
GRB comes from the activity of the central engine (per- 
haps hundreds of seconds) and the individual peaks come 
from inhomogeneities or instabilities in that activity. 

2. Expected Evolution from a Relativistic Shell 

In Fenimore, Madras, and Nayakshin (1996), we calcu- 
lated the time history and spectral evolution expected 
from a relativistic shell under the assumption of "local 
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spherical symmetry" (LSS). By local, we mean that only 

regions lying within angles the order of contribute. 
By spherical, we emphasize that the curvature of the 
shell within is the dominant feature that produces 
the envelope. By symmetry, we emphasize that the shell 
is seen head-on, rather than edge-on, as in many scenar- 
ios involving jets. Each of these assumptions is justified 
in Fenimore, Madras, and Nayakshin (1996). 

We denote the rest frame time of the detector as t 
and the arrival time as T. While the shell is gamma-ray 
active, the emitting surface keeps up with the emitted 
gamma-rays so they arrive at the detector over a much 
shorter period of time than they were emitted. If the 
material directly on axis emits for time t, the photons 
arrive on a time scale T = (c — v)t/c w t/{2r^), where 
V is the bulk velocity, T = (1 - /3^)~^/^, and /3 = v/c. 
Thus, in the case of the external shock model, the energy 
is released at the central site on a time scale of less than 
a second, and the shell expands for ^ 10^ s, emitting 
for a time scale somewhat shorter than 10^ s, but the 
photons arrive within ~ 10 s (if r = 10^). 

Let P{6, (j), R) give the rate of gamma-ray production 
on the shell as a function of spherical coordinates. Indi- 
vidual bursts might have temporal structure caused by 
variations in P as a function of {0,^); however, on av- 
erage, all off-axis regions should behave like the on-axis 
regions. Since we are considering the average properties 
of bursts in this paper, we assume P{6,^,R) is inde- 
pendent of {0,(j>). It is a single, thin shell, so R = vt 
such that P{0, 0, R) = P{t). The choice of a thin shell 
is justified by the observations. If AT is the typical 
time scale for variability in a GRB with duration Tdur, 
then the thickness relative to the size of the shell is 
Ar||/ro ~ AT(2r2Tdur)~^ (see Fenimore, Madras, and 
Nayakshin 1996). 

2.1. Temporal Structure from Emission at a Single Time 
Assume the shell expands for a time {1^) in a photon 
quiet phase and then emits for a short period of time 
(i.e., P{t) = Po<5(i — te))- In terms of arrival time, the 
on-axis emission will arrive at = t^j (2r^) and the ofi^- 
axis emission will arrive later with an expected shape 



it emits (Te). The decay phase is due to photons delayed 

by the curvature. 

GRB spectra can often be fit by the so-called "Band" 
model (Band et al. 1993) which consists of two power 
laws (indexes oix^a^) and a connecting energy. If 
E'^ is the connecting energy in the rest frame of the 
shell, the observed connecting energy is boosted, = 
Bs{T,Te)E'^. The delayed photons are boosted less than 
at the peak because they originate from regions moving 
at the angle 6 relative to the on-axis regions (see Katz, 
1994). The boost as a function of arrival time is (see 
eq. 5 in Fenimore, Madras, Nayakshin 1996) 

Bs{T,Te) = A-i = r(i -pcose) = 2r(r/re)-i (2) 

This is the most important equation in this paper since 
it establishes a one-to-one relationship between the angle 
at which photons originate and the time at which they 
arrive at the detector if Tg can be determined. It is 
from this equation that we determine how the observed 
spectra should evolve relative to how the time history 
evolves. 

2.2. Temporal Structure from Emission from a Range of 

Times 

More complex envelopes can be found from weighted 
sums of Vs{T,Te). For example, the shell might emit 
for a range of times starting at to- We model P{Te,To) 
to be non-zero from to to tmax and, we assume P{Te, To) 
can be approximated as 



P{Te,To)=Po{T,/To)-^. 



(3) 



In terms of arrival time, the on-axis emission will arrive 
between Tq = to/{2T'^) and T^ax = tmax/(2r2). The 
expected envelope, V{T), is (see eq. 11 in Fenimore, 
Madras, & Nayakshin 1996): 



v{T)= r vs{T,n)p{n,To)dn 



(4) 



Due to the curvature of the shell, off-axis photons will 
be delayed, and most emission will arrive later: 



Vs{T,T,)=0 



if T < Te (la) 



ifr>re (16) 



where wc have assumed that the rest frame photon- 
number spectrum is isotropic with a power law with in- 
dex —a, and V'l is a constant. Note in Fenimore, Madras, 
& Nayakshin, this was incorrectly dervied as a power law 
with —a — 2. This envelope is similar to a "FRED" (fast 
rise, exponential decay) where the slow, power law decay 
depends only on the time that the shell expands before 



V{T) = 



where 



if T < To (5a) 
if To < T < T„,ax (56) 

if T>Tmax (5c) 



w = a + 3 — T] 



(6) 



Here, as a simplification to show the nature of the en- 
velope, we have used a power law spectral shape with 



a number index of a (~ 1.5). The deviations from the 

Band model are not important and a ^ (ai + Q!2)/2. 

The envelope in equation (5) is also similar to a 
"FRED" where the rise depends mostly on the duration 
of the photon active phase (Tmax — ?o) and the slow, 
power law decay depends mostly on the final overall size 
of the shell (T^ax). 

The spectrum of the shell might change while it is 
emitting. We model the possible change in E'^ as 



(7) 



The observed connecting energy as a function of arrival 
time is 



Ec{T) = 
such that 



J^^ Vs{T, T,)PiT„ n)E[{T,,, To)Bs{T, T,) dT, 



V{T) 



(8) 



Eo{T) = ^ ^T^""" if To < T < T_ 



(9a) 

_ 2TE'q uj T^^^ ~ 1 if T>T (m 

-io ^ -'max -'O 

where 

u;' = u!-l'+l (10) 

3. Comparison to Observations 

In the external shock models, there is a single relativistic 
shell which should provide the overall GRB emission en- 
velope. From the temporal structure of GRBs, it is clear 
that some process is required to break the local spheri- 
cal symmetry and produce a chaotic time history rather 
than a smooth profile as given by equation (5). How- 
ever, on average, GRB time histories should conform to 
the envelope if a single relativistic shell is responsible. 
Figure la is the average time history of 35 bright long 
GRBs (from Fenimore, 1997b). It was generated by scal- 
ing the measured duration to a standard length which we 
call Tioo and then adding the 35 bursts together. (Tioo 
is an average of T50/O.5 and T90/O.9) Adding together 
in this fashion is very reasonable for the external shock 
models. Each burst appears to have a different dura- 
tion in arrival time because arrival time is compressed 
by r~^. Normalizing each burst to a standard duration 
is equivalent to normalizing to a standard F. 

The average profile is very similar to a FRED, which 
is remarkable because only 3 or 4 of the 35 bursts are 
clearly FRED-like. It is difhcult to detect the FRED- 
like envelope in individual bursts because the shell must 
be patchy in order to make a long complex burst (Fen- 
imore, Madras, & Nayakshin 1996). However, the en- 
velope becomes evident when many bursts are averaged 
together. 



The dotted line in Figure la is the best fit of equation 

(5) to the average profile. Besides an overall scale factor 
and a, there are only two free parameters: Tq and Tmax- 
Note we have assumed that F is constant. From the fit 
in Figure la, Tq = 0.78Tioo and T^ax = 0.85Tioo. We 
have shifted the T scale such that zero is at the time 
(relative to Tioo) when the central explosion occurred. 
Thus, if a relativistic shell is responsible for the shape in 
Figure la, we learn that the gamma-ray quiet phase for 
the shell is ~ l.SF^Tioo and the gamma-ray active phase 
is short: O.OTF^Tioq. 

Figure la implies that a relativistic shell can explain 
the average envelope of emission of GRBs. However, the 
single relativistic shell paradigm also predicts how the 
spectra should evolve. The dotted curve in Figure lb 
uses the fit parameters from Figure la in equation (9). 
Since there are no additional free parameters, the spec- 
tral evolution is entirely determined by the fit to the time 
history. We have averaged the 16 channel BATSE MER 
data in the same manner as the averaging in Figure la. 
The MER data were then summed into five samples of 
duration 0.15Tioo starting at T = I.OTioq. The result- 
ing spectra were fit with the "Band" spectral shape. We 
determined Ec in a manner similar to how Liang and 
Kargatis (1996) determined spectral evolution in GRB 
pulses. We first fit the entire time history to determine 
an average ai and a2- The five spectra were analyzed 
with these average power law indexes such that Ec was 
the only free parameter. In Figure lb, we have plotted 
the ratio of Ec for each sample to Ec of the first sample. 
The first sample was chosen to be in the decay phase of 
the envelope because the spectral evolution (eq. [9]) is 
effectively independent of the intrinsic variation during 
the decay phase. The average GRB envelope has spectra 
that evolve somwhat faster than a relativistic shell. 

It could be argued that the averaging of the 35 bursts 
is uncertain because each burst is chaotic, and it is dif- 
ficult to properly identify how to scale each one. Since 
FRED-like bursts have clear envelopes and are the shape 
expected from a relativistic shell, they should be easy to 
scale properly. In Figure Ic we averaged six FRED-like 
BATSE bursts, stretched by their durations to a stan- 
dard time. The decay phase fits better, probably because 
the chaotic nature of the long complex bursts only allows 
a rough determination of how much stretching is needed 
before averaging. For the FRED-like bursts, the time 
of expansion from the central site is somewhat smaller, 
To ^ O.53T100, and the rise is very rapid, implying that 
P(Te,To) is non-zero only for a short range of times. In 
Figure Id we show the expected evolution of the con- 
necting energy (dash curve) and the observed evolution 
(squares) . The conflict between the observations and the 
predictions from a single relativistic shell is even stronger 
than observed with the 35 long bright bursts. 




Fig. 1. Comparisons of time evolution expected from a single relativistic shell to GRB observations. In the "external" shock model, the overall 
envelope of emission is related to the properties of a single relativistic shell. 

(a) Fits of a single relativistic shell to the average envelope of 35 bright, long BATSE bursts. The temporal scale of each individual burst 
was scaled to a standard duration called Tioo. The dash line is the best fit solution based on equation (5). The zero on the time scale 
represents when the central explosion occurred and it is set by the fit to equation (5). The shell expanded in a gamma-ray quiet phase 
that lasted ~ 0.78Tioo and then produced gamma-rays until ~ 0.85Tioo. 

(6) Fit to the spectral evolution of the 35 bursts based on the BATSE MER data. The dash line is the expected relative characteristic 
energy in the spectrum based on the best fit parameters from (a) and equation (9). The solid squares are the observations. 
(c)Fits of a single relativistic shell to the average envelope of 6 FRED-like BATSE bursts. 

{d) The dash line is the expected evolution based on the fit in (c) and equation (9). The squares are the observed spectral evolution based 
on the BATSE MER data for the 6 FRED-like bursts. 

The spectra of GRBs evolve faster than allowed by the single relativistic shell that is utilized in the "external" shock models. The fits used 
in this figure are not substantially affected by intrinsic variation in the shell's intensity or spectra. However, variations in the bulk Lorentz 
factor might allow for a consistent fit of time history and spectra. 



4. Intrinsic Spectral Variations 

In Figure 1, we assumed no intrinsic variations in the 
emission as a function of time, that is, = 77 = 
such that P{T) = Pq and E'^ = Eg. In this section 
we show that intrinsic variations cannot provide soften- 
ing as rapid as seen in Figure 1. In Figure 2a, we use a 
range of parameters in equation (5) to demonstrate the 
range of possibilities. (All curves have Tq = 80.) The 
solid curve uses an intermediate value of Tmax (=145), 
and a value of rj of about 9 was used to obtain a fast rise 
and decay phase commensurate with the observed char- 
acteristics of GRBs. The dot-short-dash curve was found 
as in Figure 1 by fitting to the simulated data from the 
solid curve with r] = 0. The resulting Tmax is 98. The 
short-dash curve used T,„ax = 101 and rj = 5. Finally, 
the long-dash curve uses Tmax = oo and rj = 12. In all 
cases, the decay phase is approximately independent of 
the parameters used and scales roughly as (T/To)~"~^. 



Figure 2b shows the expected variation in Ec based 
on the parameters used for the curves in Figure 2a. For 
the solid curve and the dot-short-dash curve, we used no 
intrinsic variation in E'^ (i.e., v = 0). For the short-dash 
and long-dash curves, we used u = —2 and +2, respec- 
tively. During the regions where P(Te,To) is non-zero, 
the observed Ec can be influenced by intrinsic variation 
in E'^. However, in the decay phase, Ec is almost inde- 
pendent of v and varies as (T/Tq)^^ (eq. 9, Fig. 2b), 
whereas the observations (Fig. lb, d) vary faster. 

One assumption that we have implicitly made is that 

a weighted sum of Band spectra gives the appearance 
of having a value of Ec that is similar to the average of 
the values of individual Ec. We tested this assumption 
by producing simulated spectra using the parameters of 
the long-dash curve in Figure 2. Simulated spectra were 
produced for each time bin. We processed the simulated 
spectra in the same manner as the MER data and the 
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Fig. 2. Variations in the observed time history and spectra of a relativistic shell due to intrinsic variations in the rest frame of the shell. 

(a) Time histories from a range of parameters in equation (5). The key parameters are the range of times the shell emits in its rest frame 
(T() to Tmax) and how the emission evolves with time (oc T~^). The solid, dot-short-dash, short-dash, and long-dash curves correspond 
to (ro,Tmax,»;) = (80, 145, 9), (80, 100, 0), (80, 100, 5), and (80, oo, 12), respectively The intrinsic variation mostly affects the rise 
portion of the envelope. The decay phase is due to off-axis emission which arrives later due to the curvature of the shell. As a result, it is 
an average over the intrinsic variations and follows a universal shape that is due to the curvature (ex T^°^^). 

(6) Spectral evolution due to the intrinsic variation in the shell. The parameters of the curves in (a) were used in equation (9) to give 
the evolution of a characteristic energy in the spectrum. Here, Ec is the energy that connects the two power laws in the Band model. 
An additional parameter gives the intrinsic variation of Ec in the rest frame {E'^ oc T~"). The four curves used v = 0, 0, -2, and -|-2, 
respectively. We fit the Band model to simulated MER data using the parameters of the long-dash curve. The resulting Ec's are the 
squares. Thus, if one averages a series of Band spectra, the Ec of the resulting spectrum is similar to the average of the individual Ec's; 
there are no systematic effects due to the use of the Band model. 

The intrinsic variation mostly affects the spectrum in the rise portion of the envelope. The decay phase has spectra that are shifted, 
averages of the spectrum during the rise. The shift is the Lorentz factor corresponding to the location of the off-axis emission responsible 
for the decay portion. The Lorentz factor varies as , so the spectral evolution is a universal function (oc T~^) that is independent of 
the intrinsic variation. 



results are shown as the squares in Figure 2. In the decay 
phase, they closely follow the expected relationship. 

5. Discussion 

The only approximations that wc have made are that 
P and -Eg are approximately power laws (eqs. [3] and 
[7]). This allowed analytic solutions which show the na- 
ture of equations (4) and (8). These equations use only 
kinematic considerations to demonstrate that the decay 
phase is influenced only by the curvature of the shell. 
During the decay phase, the delayed ofF-axis photons 
are strongly beamed, and they are Lorentz-shifted out 
of the bandpass of the instrument. The combination of 
these two effects causes the time histories to decay as 



~ A-""i ~ (T/ro)"""^ In contrast, the characteris- 
tic energy of the spectra varies as ^ ^ (T/Tq)^^. 
For example, by the time T = 2T(), the time history will 
have returned almost to its baseline, while the charac- 
teristic energy will have changed only by a factor of 2. 
These relative rates of change are inconsistent with the 
observed overall envelope of GRB evolution. Indeed, it is 
well known that the characteristic energy of GRB spec- 
tra can change by factors the order of 4 throughout a 
burst. 

The fundamental reason that the spectra evolve as 
can be seen in equations (8) and (9). Notice that 
the only dependency on T in equation (8) comes from 
Bs{T,Te) and Vs{T,Te). To first order, the dependency 
on Vs{T, Te) is cancelled by the V{T) in the denominator 



so that EciT) - Bs{T,T^,^). If T > T^ax, then equa- 
tion (9b) states that Ec{T) — kT^^ where fc is a con- 
stant. In fact, when T > Tmax, Ec{T) is precisely a con- 
stant times for aU functions P(T, Tg) and E^iT, T^), 
regardless of their functional form. Even if P(T, Te) has 
a form that allows Tmax to formally go to infinity, there 
must be an effective T^ax where P{T, Te) approaches 
zero. Otherwise, V{T) would diverge or GRB profiles 
would not have fast rises. Thus, in all cases, Ec{T) fol- 
lows T^^ in the decay phase. 

From a physical perspective, this universal spectral 
evolution can be understood on the basis of symmetry. 
We are not claiming that the intrinsic variations in P and 
Ec are not manifested in the observations. Indeed, while 
P and Ec are changing, the observations track P and Ec 
(between T = 80 and T ~ 100 in Fig. 2). The universal 
spectral evolution only occurs in the decay phase. Dur- 
ing the decay phase, photons are delayed by the curva- 
ture of the shell, so they come from some off-axis annulus 
at an angle 6 = 9i. By symmetry, the same averaging 
that occurs at ~ also occurs at 6* = 0i; therefore, 
the two spectra are identical, except that the average 
spectrum from 6 — 9i is Lorentz shifted by the boost 
r(l — /3 cos 6*1) rather than r(l — /3). Equation (2) gives 
a one-to-one relationship of the boost to the arrival time. 
It is that boost that gives the universal spectral evolu- 
tion. 

The disagreement between the time history evolution 
and the spectral evolution is further evidence that some- 
thing must break the local spherical symmetry if single 
relativistic shells are involved in GRBs. For example, if 
the shell can act as a parallel slab or if the GRB produc- 
ing regions are confined to narrow pencil beams, then 
the envelope and spectra can evolve uncoupled. In Feni- 
more, Madras, & Nayakshin (1996), we presented several 
other arguments that implied that LSS must be broken 
or GRBs must be central engines. For example, precur- 
sors and the fact that many long complex bursts do not 
follow the expected envelope imply that many parts of 
the shell never emit photons, resulting in a low efficiency 
for converting the kinetic energy of the shell into gamma- 
rays. Sari & Piran (1997) argue that the efficiency is 
too low (~ 1%) to be consistent with merging neutron 
stars and, thus, one must resort to internal shock models 
(however, see Fenimore 1997a). 

In this paper, we show the need to break the LSS is also 
supported by the spectral evolution. We have shown that 
intrinsic variations alone are not sufficient to explain the 
observed spectral evolution. This reinforces the "shell 
symmetry" problem we introduced in Fenimore, Madras, 
and Nayakshin: GRB models that involve a single rela- 
tivistic shell must explain how the material is confined 
to pencil beams narrower than F^^, how LSS is broken, 
or how the shell can be very inefficient without excessive 



energy requirements. 

In this paper, we did not vary F with time. Varying F 
would not introduce any additional intrinsic spectral evo- 
lution since it would be equivalent to varying E'c which 
is already accounted for in equation (7). However, in a 
future paper we will show how changes in F in time can 
affect the shape of the envelope and, perhaps, result in 
the breaking of LSS. 

If the "shell symmetry" problem cannot be solved for 
the external shock models, those models might have to 
be abandoned in favor of the internal shock models, al- 
though these might have similar problems. In the in- 
ternal shock models, peaks in the GRB time history 
arise from variations in the central activity. We have 
fit equations (5) and (9) to individual peaks in GRBs 
and the spectra evolve faster than (T/Tq)^^. In gen- 
eral, characteristic energies in peaks can evolve by fac- 
tors of 4 whereas spectra from relativistic shells should 
only evolve by a factor 2 across a peak. 

If the internal shocks occur because fast shells run into 
slow shells (Rees and Meszaros 1994) resulting in a shell 
with curvature, then the internal shock models will also 
have to solve the "shell symmetry" problem. The inter- 
nal shock models also have to explain how a source can 
produce (il/47r)10^^ erg sec ~^ of activity for hundreds 
of seconds. (Here, fl/Air is the fraction of the sky into 
which the burst beams its shells.) 
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